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I. INTRODUCTION 


Ballistic projectiles intended for employment against low flying aircraft 
are characterized by a very short time of flight, which implies a flat trajectory 
and maximum target penetration/capture coupling. Current gun technology permits 
typical launch velocities to M = 5 and anticipated target defense appears 
impervious to penetration below M = 2. The long rod penetrator has limited 
applicability in this instance, being expected to pass through the target 
completely. Thus the projectile of interest is spin stabilized and operates 
at 2<M<5; and the challenge of design lies in our ability to predict the aero- 
dynamic behavior of classical shell in this environment. 


Various analytic and empirical determinations of zenon vague characteristics 
of spin stabilized projectiles are found in open literature!-° where selected 
Characteristics of the overall flight performance are considered in detail. 

The predictions of drag coefficient tend to be substantiated by test results 
while the estimates of static moment, reported less frequently, are usually 
found more difficult to verify. Generalized techniques to predict the damping 
moment and Magnus moment coefficients are seldom explored although extensive 
range data has been accumulated, 


ed McCoy* has presented an authoritative treatise on calculating drag 
coefficient and R.H. Whyte@ has published reliable data on the static moment 
coefficient. Historically, H.P. Hitchcock? has proposed a set of correlations 
of test results while Reference 4 provides an extensive bibliography with "how 
to" examples. A current analysis of the flow field around the typical projectile 
1s demonstrated by W.B. Sturek, et. al.,°»© by means of a three~dimensional 
finite difference algorithm where the damping and Magnus contributions are 
explicitly determined. These programs are available in the FORTRAN programming. 


Dee ii, MeCoy, "MeDrag" - A Computer Program for Estimating the Drag Coeffietents 
of Projectiles," BRL Report No. ARBRL-TR-02293, February 1981 (AD #A098110). 


“2, Whyte, "Spinner" - A Computer Program for Predicting the Aerodynamic 


Coeffiretents of Spin Stabilized Projectiles," General Electrie Technical 
Report TIS69APB3, 1969. 


"a P, Huteheoek, "Aerodynamte Data for Spinning Projectiles," Balltstte Research 
Laboratories Report No. 620, October 1947. (AD #800469). 


4AMCP 706-242, "Design for Control of Projectile Flight Characteristics," 1966. 


eed Sturek, D.C. Mylin and C.C. Bush, "Computational Parametric Study of the 


Aerodynamics of Spinning Slender Bodies at Supersonic Speeds," AIAA-80-1585-CP, 
August 1980. 


"ey Be Sturek, D. C. Mylin, "Computational Parametrie Study of the Magnus Effect 


on Boattatled Shell at Supersonic Speeds, AIAA-81-1900, August 1981. 


The present report offers a simplified procedure to estimate values of 
the initial aerodynamic design parameters required for flat fire trajectory 
specification within 2<M<5. Execution is by means of a desk-top computer 
program for application into the Tektronix 4051 computer and is suited for 
conversion to any comparable machine. It requires only the geometry and physical 
properties of the projectile and presents the dimensional outline, a complete 
tabulation of the aerodynamic coefficients and a graphical exposition of the 
summarized results. The aerodynamic basis is wholly empirical. 


II. PROCEDURE 


Reference 7 contains a collection of range data on a wide variety of 
projectiles which were fired in the BRL Aerodynamics Range. The maximum 
diameter of the projectiles extends to 20 mm and the geometry includes nose 
shapes of circular arc and conical development, variations in length of nose 
and afterbody, boat tail proportion and rotating bands. Figure 1 shows a line 
drawing of a typical projectile. Associated performance curves’ allow an 
insight into the effects of these variables through the supersonic flow region 
and permit the formulation of the equations proposed by this report as 
presented in Tables 1, 2, 3 and 4. One exception is that normal force and 
static moment data for conical nosed, boattailed projectiles are not available 
although the square based data is included. The performance curves are not 
fitted statistically since the survey data do not include discrete numerical | 
information nor do they indicate error bounds. However, the resulting algebraic 
expressions appear to be consistant and reasonable within the indicated Mach 
excursion. Table 5 lists the equations of the gyroscopic stability factor, 
the initial yaw period® and the velocity decrement’. The normalized caliber 
notation, where a representative linear dimension provides a reference length 
and the mass/force dimension is converted to a specific gravity, is employed 
throughout, except that the trajectory range is given metric expression. 


A complete program listing with initialization instructions is presented 
in Appendix A. By command print-out, the program provides 


e An outline plot of the projectile. 
e A complete tabulation of aerodynamic coefficients. 


e Individual graphs of the drag, normal force, static moment, damping 
moment and Magnus moment slopes versus Mach number. 


e A graph of velocity versus range with annotated schedule of retardation. 


e The length of the initial yaw cycle and a determination of the gyros- 
ecopic stability f£acter. 


WB Braun, "Aerodynamte Data for Small Arms Projectiles," Ballisties 


Laboratortes Report 1630, January 1973, (AD #909757L). 


SCLH. Murphy, "Free Flight Motion of Symmetrte Missiles," BRL Report No. 1216, 
July 1963 (AD #442757). 


Wi Donovan, "Simplified Determination of Retardation for Kinettie Energy 


Projeettles," BRL Memorandum Report No. ARBRL-MR-02994, February 1980 (AD # 
086095). 
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ITI. RESULTS 

With test data to M = 3.2, the HS831-L projectile as described by Neitzel?? 
is used to demonstrate the working procedure. The computer program (Appendix 
A) lists the input, Table A, and also develops the output according to the 
following schedule. 


Figure No. Description 

2 Projectile Outline 
(Table 6) Table of Aerodynamic 

. Coefficients vs M 

5 : Cyr VS M 

4 Cy VSM 

5 orn vs M 

6 Cig Cus vs M 

7 OM o vs M 

8 Trajectory Parameters 


Test values from Reference 10 are superposed on the graphs by asterisks. The 
predictive code indicates conservative values for this projectile. 


IV. CONCLUSIONS 


With the facility of the developed programs as employed by the desk-=top 
computer, the design problem of estimating values of the aerodynamic coefficients 
and trajectory parameters for spin-stabilized shell in flat fire is considerably 
simplified. Individual effects of geometric modification of a particular 
candidate projectile can be determined immediately. Expertise in the science 
of fluid mechanics is not required to produce rational results, simple 
constraining criteria being sufficient for initial considerations. 


Although a number of projectiles has been examined by these programs, 
only the M-392!1 and those of Reference 5 show data beyond M = 3.6. The M-392 
1s unusual in that a base stub is attached (Appendix B) while the Reference 5 
models, Appendix C) are precision manufactured for wind tunnel study. Neither 
carries a rotating band. 


tO ge, Nettzel, Jr., "Aerodynamic Characteristies of 30 MM HS831-L Ammunition 
Used in the British 30 M Rarden Gun," Ballistie Research Laboratories 


Memorandum Report 2466, March 1975 (AD #B003939L) 


11 My. Sehmidt, B.P. Burns and G. Samos, "Replica Modeling of the Launch and 


Flight Aerodynamtes of Projectiles," BRL Report No. ARBRL-TR-02104, 
September 1978 (AD # A063521). 


17 


The M-392 determinations are included, although these represent a unique 
projectile, to indicate the very strong effect of the base stub on the normal 
force and static moment development. If this constitutes a general behavior 
pattern, and such a conclusion requires correlating verification, the phenome- 
nology must be examined by analytical techniques. 
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AERODYNAMIC COEFFICIENTS FOR SPIN STABILIZED PROJECTILES 


TABLE 6. 
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APPENDIX A 
INITIALIZATION INSTRUCTIONS 


Zo 


(1) The program is stored using three files of a magnetic tape cartridge 
intended for use on a Tektronix 4051 desk-top minicomputer. To begin program 
execution, press the "'autoload" button or use the dual commands of ''Find 1" 
followed by '"'Old''. It will take approximately 15 seconds to load the first 
program and from this point on, the user is guided through the program 

by a number of questions. 


(2) The program consists of two main subprogram files, one for the conical 
nose shape and the other for the ogival nose shape. The user is asked the 
desired shape and is taken to the respective file. Both subprograms offer 
identical menus but involve different calculations. 


(3) A maximum of 20 physical parameters is required for program execution. 
All such parameters must be input in calibers with a few noted exceptions. 
These parameters are requested as follows: 


a) nose length 
b) length of radius 
Cc) cylindrical body length 
d) boattail length 
e) boattail angle 
= if} radius of curvature offset 
g) rotating band radius 
h) rotating band width 
1) distance from rotating band center to base of projectile 
AD length of boom extension 


*--QOnly required for ogival nose 


The graph of velocity versus range requires eight additional input 
parameters. Because this graph is not always desired, receiving it is 
optional to the user. The additional parameters are requested as follows: 


a) axial moment of inertia 

b) transverse moment of inertia 

c) weight of projectile 

d) twist of projectile 

e) diameter of projectile in millimeters 
£) air density 

g) muzzle Mach number 

h) maximum range of projectile desired. 


Both subprograms use thirteen User Definable Keys. These thirteen selections 
are listed below. 


. Nomenclature 
..Table of Aerodynamic Coefficients 
..Total Drag Coefficient (C, ) versus Mach Number Plot 


..Total Normal Force Slope cSefficient CC 2 ) versus Mach Number Plot 
..Static Moment Slope Coefficient Coy versus Mach Number Plot 
..Damping Derivative Slope Coefficient (C,,-+C,, ) versus Mach Number Plot 


..Magnus Moment Slope Coefficient (Choa Versi Mach Number Plot 


oO NDS WN 


.. Velocity versus Range Plot, Gyroscopic Stability Factor, 
Initial Yaw Cycle, Retardation at different range values 


alt 


..Return to original spin stabilized projectile design 
.. Views the initial data (see (5)) 

..Calculates data (see (5)) 

..Prepares the user to view another projectile shape 
--Returns user to menu 
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“,. DESIGN | 
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Flowchart for Automatic Plotting Routines 
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Intercept of C 


Slope of C 


D 


LIST OF SYMBOLS 


n vs M characteristic 


vs M characteristic 


Center of gravity of projectile 


Center of pressure of normal force 


1.0 caliber, reference diameter 


Afterbody length 


Boattail length 


Nose length 


Length of rotating band moment arm 


Length of base extension 


Mass of projectile 


Rotation rate 


Nondimensional spin rate 


Transverse angular velocities 


(o 


Z 


Ll. 
eG 


Radius of face of base 


Range 


Velocity of projectile 


Sonic velocity of ambient atmosphere 


Distance from nose to center of gravity 


Distance from nose to center of pressure of normal force 


Reference area, (.785 cal) 


4pv A 


Pressure drag coefficient, - base of body 


Z 


Drag Force 


ref 


Rotating band drag coefficient 
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“Na (bt) 


Na(ext ) 


I 
x 


LIST OF SYMBOLS (continued) 
Viscous drag coefficient - body 
Wave drag coefficient - body (nose) 


Total drag coefficient 


Static Moment 


2 
1. 


,» aerodynamic moment slope coefficient 


Aerodynamic moment slope with base extension 


Magnus Moment , Magnus moment slope coefficient; moment 


2 d e e »° 
1. pe 
3 pv A nef ( eo 6 rotates nose | to plane of O in direction 
of spin. 
ramping “onent , damping moment slope coefficient; moment 
4% OV A nef d(4t ) increases angular velocity 
oe se » normal force slope coefficient 
1 
4p Vv A 6 
ref 


Normal force slope coefficient for projectiles with boattails. 
Normal force slope coefficient for projectiles with base extensions. 
Axial mass moment of inertia 

Transverse mass moment of inertia 


Mach number 


Mach number at muzzle 


cM. + b 
0 
Operational parameter, M 
0 
ref 2 
Operational parameter, = 


Ogive nose radius 


Angle of attack, sideslip (employed as a directional subscript) 
L. 
, total angle of attack 


Rate of change of a 


Air density 


Sin 
eo, 
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Please take a few minutes to answer the questions below; tear out 
this sheet, fold as indicated, Staple om tape closed, and place 
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improving future reports. 


1. BRL Report Number 
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